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Key Points: 


¢ The role of solar forcing in global warning is unclear due to the uncertainty in the 
change in TSI since the Maunder minimum 

¢ The TSI level of the Sun when it is in its least active state is established from an 
advanced model of TSI variability 

¢ This lower limit on grand minima TSI truncates the current range of the possi- 
ble change in TSI since the Maunder minimum by more than half 
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Abstract 

How the solar electromagnetic energy entering the Earth’s atmosphere varied since pre- 
industrial times is an important consideration in the climate change debate. Detrimen- 
tal to this debate, estimates of the change in total solar irradiance (TSI) since the Maun- 
der minimum, an extended period of weak solar activity preceding the industrial revo- 
lution, differ markedly, ranging from a drop of 0.75 W m7? to a rise of 6.3 W m7”. Con- 
sequently, the exact contribution by solar forcing to the rise in global temperatures over 
the past centuries remains inconclusive. Adopting a novel approach based on state-of- 
the-art solar imagery and numerical simulations, we establish the TSI level of the Sun 
when it is in its least active state to be 2.00.7 W m~? below the 2019 level. This means 
TSI could not have risen since the Maunder minimum by more than this amount, thus 
restricting the possible role of solar forcing in global warming. 


Plain Language Summary 


How the amount of energy the Earth receives from the Sun varied since pre-industrial 
times is an important consideration in the climate change debate. Detrimental to this 
debate, it is not known whether the Sun grew brighter or dimmer since the 16th century, 
and by how much. As a consequence, the exact contribution by fluctuations in the bright- 
ness of the Sun to the rise in global temperatures over the past centuries remains con- 
troversial. It is established that the Sun was particularly inactive over much of the 16th 
century. Adopting a novel approach based on state-of-the-art solar imagery and com- 
puter models, we determined the brightness of the Sun when it is in its least active state 
possible. This places a strict limit on how much the Sun could have grown brighter since 
the lull in solar activity over the 16th century, restricting the possible role the Sun could 
have played in global warming. 


1 Introduction 


Climate change is one of the most pressing challenges facing humanity. The urgency 
is underscored by the recent pronouncement by the Intergovernmental Panel on Climate 
Change that greenhouse gas emissions will have to be halved by 2030 to keep global warm- 
ing below the critical threshold of 1.5 °C (IPCC, 2018). Numerical simulations of the Earth’s 
climate are central to the understanding of past climate and predictions of future trends. 
They require, as input, the time variation in the various agents of climate change, both 
anthropogenic and natural (IPCC, 2013). A key natural driver is the variability in the 
electromagnetic energy reaching the Earth from the Sun (Gray et al., 2010). This is usu- 
ally described in terms of solar irradiance, which is the solar radiative flux above the Earth’s 
atmosphere at the mean Earth-Sun distance. The long-term or secular trend in solar ir- 
radiance since pre-industrial times, and how this variability is distributed in wavelength 
are of particular interest for their relevance to the overall influence of the Sun on the Earth’s 
climate over the past centuries (Solanki et al., 2013). In this article, our focus will be 
on the particular issue that current estimates of the secular change in the wavelength- 
integrated total solar irradiance (TSI) since the 16th century diverge significantly, rang- 
ing from a drop of 0.75 W m7? to a rise of 6.3 W m~? (Tapping et al., 2007; Steinhilber 
et al., 2009; Schrijver et al., 2011; Shapiro et al., 2011; Judge et al., 2012; Coddington 
et al., 2016; Dasi-Espuig et al., 2016; Egorova et al., 2018; Wu et al., 2018; Lockwood 
& Ball, 2020). 


Since the reliable measurement of solar irradiance only started in 1978 (Kopp, 2014), 
climate simulations rely on models of solar irradiance variability to provide the requi- 
site historical solar forcing input. Except at the shortest and the longest timescales, where 
oscillations, convection and stellar evolution become relevant, solar activity is dominated 
by its magnetism, making it the prime candidate driver of solar irradiance variability at 
intermediate timescales (Solanki et al., 2013). Indeed, various studies have demonstrated 


the variation in solar irradiance at timescales of days to decades to be mainly the result 
of solar surface magnetism (Yeo, Krivova, Solanki, & Glassmeier, 2014; Yeo et al., 2017; 
Shapiro et al., 2017). This is widely considered to also be the case at longer timescales 

of centuries and millennia (Solanki et al., 2013). The photospheric magnetic field is par- 
tially confined in kilogauss-strength magnetic concentrations, which manifest themselves 
as dark sunspots and bright faculae and network (Spruit & Roberts, 1983). Solar irra- 
diance fluctuates with the changing prevalence and distribution of these bright and dark 
magnetic structures on the Earth-facing solar hemisphere. Based on this knowledge, his- 
torical solar irradiance is reconstructed from proxies of solar magnetism by modelling 

the effect of photospheric magnetism on solar irradiance (Tapping et al., 2007; Steinhilber 
et al., 2009; Shapiro et al., 2011; Coddington et al., 2016; Dasi-Espuig et al., 2016; Egorova 
et al., 2018; Wu et al., 2018). 


The various reconstructions of TSI extending back to pre-industrial times indicate 
divergent secular variability. The Maunder minimum is a solar grand minimum, a multi- 
decade period of weak solar activity, which occurred just before the industrial revolu- 
tion (around 1645 to 1715, Fig. 1A). According to the CHRONOS model (Egorova et 
al., 2018), between the Maunder minimum and the ongoing solar activity cycle minimum, 
TSI rose by 3.9 to 6.3 W m~? (Fig. 1B). This is a recent revision of the model by Shapiro 
et al. (2011), which estimated a rise of 643 W m~?. The revision was prompted by the 
assertion by Judge et al. (2012) that Shapiro et al. (2011) might have over-estimated this 
quantity by a factor of about two. The other model-based estimates (Tapping et al., 2007; 
Steinhilber et al., 2009; Coddington et al., 2016; Dasi-Espuig et al., 2016; Wu et al., 2018) 
are markedly weaker and relatively consistent with one another, indicating a rise of be- 
tween 0.3 and 0.9 W m~? (Fig. 1B). In another recent study, Lockwood and Ball (2020) 
concluded, from an analysis of the measurements and models of TSI over the past two 
solar cycles, that it cannot be excluded from the uncertainty in this data that TSI might 
have been 0.4 W m~? lower to 0.75 W m~? higher during the Maunder minimum. Fur- 
thermore, there is a claim that the last cycle minimum in 2008 might be representative 
of the Maunder minimum (Schrijver et al., 2011). This implies that TSI rose since the 
Maunder minimum by the same margin as it did since 2008, about 0.1 W m~?. By cur- 
rent estimates, since the Maunder minimum, TSI could have risen by more than 3 W 
m~? (Shapiro et al., 2011; Judge et al., 2012; Egorova et al., 2018) or by less than 1 W 
m~? (Tapping et al., 2007; Steinhilber et al., 2009; Schrijver et al., 2011; Coddington et 
al., 2016; Dasi-Espuig et al., 2016; Wu et al., 2018), or even declined instead (Lockwood 
& Ball, 2020). 


While there is clear evidence that greenhouse gas emissions play a significant role 
in global warming (Shangguan et al., 2019), quantifying its exact effects is less than straight- 
forward due to, amongst other factors, the uncertainty in the secular variability in so- 
lar forcing. For certain climate models (Feulner, 2011; Ferndndez-Donado et al., 2013; 
Hind & Moberg, 2013; Schurer et al., 2014; Luterbacher et al., 2016), assuming a 3 to 
6 W m`? rise in TSI since the Maunder minimum produces a 0.3 to 0.4 °C greater in- 
crease in Northern Hemisphere temperatures over this period than if a rise in TSI of less 
than 1 W m~? were considered. A steep secular rise in TSI since the Maunder minimum 
points to solar forcing effecting a greater temperature rise over the past centuries and 
hence a smaller role by greenhouse gas emissions in driving climate change. It also al- 
ludes to a greater offset to global warming by the Sun if it were to enter a grand min- 
imum state in the future. Unfortunately, due in part to the divergence between differ- 
ent climate models, it remains inconclusive from climate modelling if weak or strong so- 
lar forcing is more consistent with historical temperature reconstructions (Ljungqvist et 
al., 2019). There is a vital need to constrain TSI secular variability more tightly and to 
do so in a manner that circumvents the limitations of the earlier studies which led to the 
nagging uncertainty in this quantity. 
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Figure 1. (A) Solar activity since the Maunder minimum, the period within the vertical 
dashed lines, as indicated by the group sunspot number. The original time series by Hoyt and 
Schatten (1998) (black) and the recent revision by Svalgaard and Schatten (2016) (red) are de- 
picted. (B) The reconstruction of TSI over the same period from the SATIRE-T (Wu et al., 
2018) (black), NRL (Coddington et al., 2016) (blue) and CHRONOS models (Egorova et al., 
2018) (red), illustrating the uncertainty in the change in TSI since the Maunder minimum. 
(There are multiple CHRONOS TSI reconstructions, differing by the input data used. The time 
series indicating the weakest and strongest secular variability are depicted.) (C) The reconstruc- 
tion of TSI since 2010 from the current study (black) and the measurements from SORCE/TIM 
(Kopp et al., 2005) (red), which monitored TSI over the same period except between August 
2013 and February 2014. The root-mean-square difference between the two daily time series is 
0.067 W m~? and the Pearson’s correlation coefficient is 0.987, i.e., the model replicates about 
97% of the observed variability in the TIM record. We establish, with this model, that the TSI 
level of the Sun in its least active state is 2.00.7 W m~? below the 2019 level, indicated by 
the shaded range in B) and C). Since this is the lowest level TSI can take over solar grand min- 
ima, it implies TSI could not have risen since the Maunder minimum by more than this margin, 
excluding the acute secular variability indicated by CHRONOS. 


To address this exigency, we make use of the recently-developed, novel model of 
solar irradiance variability reported by Yeo et al. (2017) to constrain the rise in TSI since 
the Maunder minimum. The Yeo et al. (2017) model is advanced in that it is the only 
one reported in the literature to not require any calibration to measured solar irradiance 
variability. Instead of estimating the change in TSI since the Maunder minimum by re- 
constructing TSI back to this time, the approach and aim here is to set an upper limit 
by establishing the TSI level corresponding to the Sun in its least active state. As this 
represents the lowest level TSI can take during grand solar minima, the difference to the 
present level represents the greatest possible rise in TSI since the Maunder minimum. 


2 The Least Active State of the Sun 


Let us consider the question of what is the least active state of the Sun. Except 
at the shortest and longest timescales, solar activity is dominantly driven by its mag- 
netism. The solar magnetic field is perpetuated by global and small-scale turbulent dy- 
namo processes. The global dynamo refers to the physical processes that produce the 
cyclic regeneration of the large-scale solar magnetic field that underlies the solar cycle, 
including the manifestation of sunspots, faculae and network (Charbonneau, 2020). The 
small-scale turbulent dynamo or SSD denotes the interaction between solar convection 
and magnetic flux that yields the ubiquitous small-scale turbulent magnetic field on the 
solar surface, the so-termed internetwork magnetic field (Borrero et al., 2017). There is 
evidence to suggest that the internetwork magnetic field might also contribute to the mag- 
netic network. 


The dearth of sunspots during the Maunder minimum (Fig. 1A) implies that the 
global dynamo is weak during this and other grand solar minima. At the same time, re- 
cent studies indicate that the internetwork magnetic field is unaffected by larger-scale 
magnetic structures (Lites, 2011; Rempel, 2014) and does not vary with the solar cycle 
(Buehler et al., 2013; Lites et al., 2014), implying that the SSD is not or only weakly cou- 
pled to the global dynamo. We emphasize here that while existing studies found no clear 
evidence that the SSD might be coupled to the global dynamo, this possibility cannot 
be excluded. Under the assumption that the SSD is, as the cited studies suggest, not or 
only weakly coupled to the global dynamo, the diminishing strength of the global dy- 
namo over grand solar minima has limited effect on the SSD. The Sun evolves on timescales 
of gigayears (Sackmann et al., 1993) and the thermal relaxation timescale of the convec- 
tion zone is in the order of 10° years (Spruit, 1982). This means the convection zone, 
the engine of the SSD, is in essentially the same state today as during the Maunder min- 
imum, or indeed over the past millennia. In summary, one, the global dynamo is weak 
during grand solar minima, two, the SSD is untethered to the global dynamo and there- 
fore unaffected by any weakening of the latter, and three, solar convection, which drives 
the SSD, is invariant over the past millennia. 


All considered, the scenario where the global dynamo is dormant, leaving just the 
SSD, represents the most inactive state the Sun can be in. In this state, there would be 
no sunspots or faculae due to the dormancy of the global dynamo, and the internetwork 
magnetic field, including the network it sustains, would extend over the entire solar sur- 
face. Assuming the SSD is, as the studies we have discussed here suggest, untethered to 
the global dynamo and invariant over the past millennia, the internetwork magnetic field 
would appear as it does today. 


3 The Brightness of the Sun in its Least Active State 


To establish the TSI level of the Sun in its least active state, and therefore the max- 
imum possible secular rise since the Maunder minimum, we build on the recent, ground- 
breaking work by Yeo et al. (2017) which presented the only reported model of solar ir- 
radiance variability to not require any calibration to measured variability. This advance, 


critical to our objective here, was achieved by deriving the intensity of certain solar sur- 

face features from three-dimensional (3D) magnetohydrodynamic (MHD) simulations of 

the solar surface and atmosphere. Such simulations have matured to such a level of phys- 

ical realism that they are able to reproduce a wide range of observational constraints (Shelyag 
et al., 2004; Schüssler & Vogler, 2008; Danilovic et al., 2010, 2013, 2016; Afram et al., 

2011; Rempel, 2014; Riethmiiller et al., 2014; Yeo et al., 2017; del Pino Aleman et al., 

2018). 


Other models, for simplicity, either obtain the intensity of solar surface features from 
semi-empirical plane-parallel model solar atmospheres (Shapiro et al., 2011; Dasi-Espuig 
et al., 2016; Egorova et al., 2018; Wu et al., 2018) or estimate their effect on solar irra- 
diance by the linear regression of certain indices of solar magnetism to solar irradiance 
measurements (Tapping et al., 2007; Steinhilber et al., 2009; Coddington et al., 2016). 

As a direct consequence, they feature free parameters that have to be constrained by op- 
timising the agreement between modelled and measured solar irradiance variability. This 
renders the model output, including the reconstructed secular variability, susceptible to 
the uncertainty in solar irradiance measurements. This is compounded by the fact that 
the linear regression of indices of solar magnetism to solar irradiance measurements is 

a crude simplification of the nuanced relationship between solar surface magnetism and 
irradiance variability (Yeo, Krivova, & Solanki, 2014; Yeo & Krivova, 2019). Also, plane- 
parallel model atmospheres, the one-dimensional simplification of the spatially inhomo- 
geneous solar atmosphere, cannot capture all the relevant physics (Uitenbroek & Criscuoli, 
2011; Holzreuter & Solanki, 2013). 


In an advance, the 3D MHD simulations-based model presented by Yeo et al. (2017) 
reproduces measured TSI variability without any of these simplifications or calibration 
to the latter. The model does, however, make the assumption that the solar surface out- 
side of sunspots, faculae and network (i.e., the internetwork) is field-free, neglecting the 
SSD and the internetwork magnetic field it maintains. The results of recent simulations 
(Rempel, 2020), indicating that the brightness of the internetwork changes with the strength 
of the internetwork magnetic field, suggest that this assumption might be problematic. 
In this study, we extend the Yeo et al. (2017) model to include a physically realistic de- 
scription of the internetwork by using an updated 3D MHD code to generate simulations 
of the SSD. This not only renders the model even more physical, but critically, it enables 
us to compute the TSI level of the Sun in its least active state, where the global dynamo 
is dormant, leaving just the SSD active (Sect. 2). And the unique strength of the Yeo 
et al. (2017) approach allows us to do so while circumventing the limitations of other mod- 
els, which contributed to the gross uncertainty in TSI secular variability. 


Based on the assumption that the variation in TSI at timescales greater than a day 
is dominantly-driven by solar surface magnetism, the Yeo et al. (2017) model describes 
the variation in TSI due to the changing prevalence and distribution of bright and dark 
magnetic structures on the solar disc. The model has two inputs. The first is the solar 
disc coverage by magnetic features. These come in the form of faculae and network, classed 
together and termed collectively as faculae, and of sunspots. The disc coverage by fac- 
ulae and sunspots is derived by identifying these solar surface features in full-disc lon- 
gitudinal magnetograms and continuum intensity images from the HMI instrument on- 
board the SDO mission, launched in 2010 (Scherrer et al., 2012). The solar disc outside 
of sunspots and faculae (including network) is taken to be the internetwork. The sec- 
ond input is the emergent intensity from sunspots, faculae and the internetwork. Sunspot 
intensity is calculated from plane-parallel model atmospheres (Unruh et al., 1999), and 
facular and internetwork intensity from 3D model atmospheres using radiative transfer 
codes. Sunspots can be reasonably well-represented by plane-parallel model atmospheres 
such that their application to the model reproduces the effect of sunspots on TSI while 
maintaining its independence from calibration to measured solar irradiance variability. 
The 3D model atmospheres are based on MHD simulations of facular and internetwork 


regions executed with the MURaM code (Vogler & Schiissler, 2007; Rempel, 2014, 2020). 
Importantly, while Yeo et al. (2017) had represented the internetwork with field-free sim- 
ulations, neglecting the SSD and the internetwork magnetic field, we make use of a more 
recent version of the MURaM code than was used in this earlier study to generate sim- 
ulations that emulate SSD action and replicate the internetwork magnetic field. 


There are four simulation runs (Fig. 2A), setup to emulate SSD action (Vogler & 
Schiissler, 2007; Rempel, 2014, 2020). The first run, representing the internetwork, is ini- 
tiated with a weak, random magnetic field (root-mean-square field strength of 1073 G). 
This seed field is amplified and reorganised in the simulation into a small-scale turbu- 
lent field resembling the internetwork magnetic field. The other three runs, represent- 
ing facular regions with varying levels of magnetic flux, were initiated with a uniform 
vertical magnetic field of 100 G, 200 G and 300 G, respectively. As in the first run, these 
runs emulate SSD action and replicate the internetwork magnetic field. However, hav- 
ing imposed net magnetic flux, these runs also form kilogauss-strength magnetic aggre- 
gations in the intergranular lanes, i.e., network and faculae-like features. Ten samples 
are taken of each run and the emergent intensity from each simulation snapshot, essen- 
tially a 3D model atmosphere, is computed (Fig. 2B). Internetwork intensity is given by 
the mean intensity of the SSD-only run snapshots, and facular intensity by the inten- 
sity of the facular features apparent in the snapshots of the facular runs. MURaM sim- 
ulations are well-suited for the current purpose, having been demonstrated in multiple 
studies to reproduce the observed properties of the internetwork magnetic field (Schiissler 
& Vogler, 2008; Rempel, 2014; Danilovic et al., 2010, 2016; del Pino Aleman et al., 2018) 
and of faculae and network (Shelyag et al., 2004; Afram et al., 2011; Danilovic et al., 2013; 
Riethmiiller et al., 2014; Yeo et al., 2017). 


Given the surface coverage by sunspots and faculae at a particular time, the model 
computes the corresponding TSI level by assigning to each point on the solar disc the 
appropriate intensity depending on whether it is in a sunspot, in faculae or the internet- 
work. Using an HMI longitudinal magnetogram and the concurrent continuum inten- 
sity image from each day between 30 April 2010, when the instrument started regular 
operation, and 31 December 2019 as input into our TSI model, we reconstruct the daily 
TSI over the intervening period (Fig. 1C). The reconstruction closely reproduces mea- 
sured TSI variability. The Sun is at its least active when the global dynamo is dormant 
and the internetwork magnetic field maintained by the SSD extends over the entire so- 
lar surface (Sect. 2). And since the SSD is likely untethered to the global dynamo and 
invariant over the past millennia, the TSI level of the Sun in this state emerges from the 
model by taking the entire solar disc to resemble the present-day internetwork. The ab- 
solute scale of the model is set to that of the SORCE/TIM TSI record (Kopp et al., 2005), 
the radiometry of which is widely regarded to be the most reliable (Fehlmann et al., 2012). 
To this end, we normalise the reconstruction of daily TSI over the last decade to the TIM 
record and apply the same normalisation factor to the TSI level of the Sun in its least 
active state, arriving at the final value of 1358.740.7 W m~? (shaded range, Figs. 1B 
and 1C). We emphasize that this normalisation merely sets the absolute scale of the model 
to that of the TIM record, the TSI level of the Sun in its least active state relative to 
the present level is unaffected. The associated error, 0.7 W m~?, denotes the uncertainty 
in this quantity from the uncertainty in the simulation of the internetwork magnetic field. 
This is a conservative estimate derived by repeating the above analysis for the extreme 
scenario that the internetwork is field-free. 


Since the TSI level of the least active Sun, 1358.7+0.7 W m7? is the lowest level 
TSI can take during grand solar minima, the difference to the 2019 level, 2.00.7 W m~? 
(shaded range, Figs. 1B and 1C), represents the maximum possible rise in TSI since the 
Maunder minimum. The TSI model and the above analysis are further detailed in the 
Supporting Information. 
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Figure 2. 3D MHD simulations of the solar surface and atmosphere, and the calculated emer- 
gent intensity. (A) Magnetic flux density at optical depth unity (i.e., at the solar surface) in a 
snapshot of each simulation run and (B) the corresponding bolometric image, calculated from 
each snapshot using a radiative transfer code. The simulation runs are set-up to emulate the 
small-scale turbulent dynamo, producing the ubiquitous small-scale turbulent magnetic field. Into 
three of the four runs, a uniform vertical magnetic field of 100 G, 200 G and 300 G, respectively 
was introduced. After these simulations reached statistical equilibrium, the imposed net magnetic 
flux aggregates in the intergranular lanes, forming the kilogauss-strength magnetic concentrations 
apparent as comparatively small bright features in the bolometric images. The white box in the 
left column denotes 1x1 arcsec, corresponding to an area of approximately 725x725 km on the 


Sun. 


The occurrence of numerous and/or sizable sunspots can produce large dips in TSI, 
even below the level of the least active Sun, 1358.740.7 W m~?, such as in October 2014 
(Fig. 1C). It is unlikely, however, that sunspots might have suppressed the secular trend 
in TSI over the Maunder minimum below this level. Firstly, there was a dearth of sunspots 
during the Maunder minimum (Fig. 1A). Secondly, darkening by a sunspot is a short- 
term phenomenon and not a state the Sun can sustain over the duration of a solar cy- 
cle, much less a grand minimum. 


4 Comparison to Other Estimates of the Change in TSI Since the Maun- 
der Minimum 


The cap on the rise in TSI since the Maunder minimum derived here, 2.00.7 W 
m~?, encloses all the current estimates (Tapping et al., 2007; Steinhilber et al., 2009; Schri- 
jver et al., 2011; Shapiro et al., 2011; Judge et al., 2012; Coddington et al., 2016; Dasi- 
Espuig et al., 2016; Egorova et al., 2018; Wu et al., 2018; Lockwood & Ball, 2020) ex- 
cept the range of 3 to 6.3 W m~? indicated by Shapiro et al. (2011), the reassessment 
of their model by Judge et al. (2012), and CHRONOS (Egorova et al., 2018) (Fig. 1B). 
We consider the upper limit established here to be more reliable than the existing es- 
timates. As argued in Sect. 3, the Yeo et al. (2017) modelling approach circumvents the 
limitations of the simpler models employed in the earlier works. And as we will argue 
next, the model atmospheres employed by Shapiro et al. (2011) and CHRONOS are ill- 
suited for estimating TSI secular variability via their methodology. In these two mod- 
els, the intensity of the quiet Sun is assumed to vary with solar activity between two de- 
fined limits, and this forms the main contribution to the reconstructed secular variabil- 
ity. (The quiet Sun refers to the solar surface outside of sunspots and faculae, i.e., net- 
work and internetwork.) In Shapiro et al. (2011), the lower and upper limits are given 
by the FAL-A and FAL-C plane-parallel model solar atmospheres (Fontenla et al., 1999), 
respectively. The FAL-A and FAL-C models correspond to quiet Sun regions with dif- 
ferent amounts of magnetic flux and therefore different intensities. Judge et al. (2012) 
raised questions over the use of the FAL-A model as the lower limit, leading CHRONOS 
to employ an amalgamation of the FAL-A and FAL-C models, termed PW-B, as the lower 
limit instead (Egorova et al., 2018). 


The FAL model atmospheres are, however, not designed for their intended purpose 
in Shapiro et al. (2011) and CHRONOS. As noted earlier, plane-parallel representations 
of the spatially inhomogeneous solar atmosphere cannot capture all the relevant physics. 
Multiple studies demonstrated that plane-parallel model atmospheres cannot fully re- 
turn the true average intensity spectra of the regions they represent (Uitenbroek & Criscuoli, 
2011; Holzreuter & Solanki, 2013). More critically, while solar electromagnetic radiation 
and consequently TSI variability is produced mainly in the photosphere, the photospheric 
layers of the FAL model atmospheres, and therefore also those of PW-B, are poorly con- 
strained (Vernazza et al., 1981). This is due to the fact that the FAL model atmospheres 
are based on ultraviolet measurements formed largely in the chromosphere. As an ex- 
ample, if the PW-B model were to be made a mere 10 K warmer, a margin smaller than 
the accuracy of this and other such model atmospheres, CHRONOS TSI secular vari- 
ability would decrease by a factor of about five (see Supporting Information). Consequently, 
the FAL models are not accurate enough to support the acute TSI secular variability in- 
dicated by Shapiro et al. (2011) and CHRONOS. The approach taken here, by relying 
on 3D model atmospheres based on realistic simulations of the solar surface and atmo- 
sphere that reproduces a wide range of photospheric observational constraints (Shelyag 
et al., 2004; Schüssler & Vogler, 2008; Danilovic et al., 2010, 2013, 2016; Afram et al., 
2011; Rempel, 2014; Riethmiiller et al., 2014; Yeo et al., 2017; del Pino Alemán et al., 
2018), circumvents these limitations. 


5 Conclusion 


There is significant uncertainty in the historical secular variability in solar irradi- 
ance, with current estimates of the change in TSI since the Maunder minimum ranging 
from a drop of 0.75 W m~? to a rise of 6.3 W m~?. This in turn introduces consider- 
able uncertainty in the contribution of the Sun to climate change. By demonstrating that 
the Sun in its least active state is at the most 2.0+0.7 W m~? dimmer than it was over 
2019, we greatly constrain the change in TSI since the Maunder minimum by ruling out 
rises above this margin. This constraint on TSI secular variability over the past centuries, 
by limiting the possible contribution by solar forcing to global temperature changes, will 
be a boon to the climate change assessment. 


Acronyms 


3D 3-Dimensional 

CHRONOS Code for the High spectral ResolutiOn recoNstructiOn of Solar irradiance 

FAL-A Fontenla, Avrett and Loeser model A 

FAL-C Fontenla, Avrett and Loeser model C 

HMI Helioseismic and Magnetic Imager 

IPCC Intergovernmental Panel on Climate Change 

MHD MagnetoHydroDynamic 

MURaM Max Planck Institute for Solar System Research and University of Chicago 
Radiation MHD code 

NRL Naval Research Laboratory 

PW-B Present Work model B 

SATIRE-T Spectral And Total Irradiance REconstruction for the Telescope era 

SDO Solar Dynamics Observatory 

SORCE SOlar Radiation and Climate Experiment 

SSD Small-Scale Dynamo 

TIM Total Irradiance Monitor 

TSI Total Solar Irradiance 
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